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Extensive ARPES and low-energy inelastic neutron scattering studies of cuprate superconductors 
can be successfully described using a weak-coupling theory in which quasiparticles on a square lattice 
interact via scalar and spin-dependent effective interactions. In this article we point out that in un- 
derdoped Bi2212 both probes are consistent with dominant near-neighbour Heisenberg interactions. 
We discuss the implications of this finding for the mechanism of high- Tc superconductivity. 

PACS numbers: 



I. INTRODUCTION 

It is evident from many experiments^i^i^ that the su- 
perconductivity of the high- Tc cuprates can be described 
at low energies and temperatures by a BCS theory with 
effective interactions between square-lattice quasiparti- 
cles that lead to short-coherence-length d-wave supercon- 
ductivity. After many years of study, the source of this 
effective interaction has still not been established with 
certainty. The d-wave property is naturally associated 
with near-neighbor interactions, but these could be spin- 
independent and attractive (V) or Heisenberg-like and 
antiferromagnetic (J). The effective interactions could 
be mediated by a lattice or electronic-fluctuation boson, 
as in conventional superconductivity, or fall outside of the 
familiar Eliashberg scenario. The resonance feature in in- 
elastic neutron scattering^, which appears to be generic 
in cuprates but absent in conventional superconductors, 
can be explained^i^iii^ if interaction parameters are cho- 
sen so that the system is close to an antiferromagnetic 
instability, possibly one driven by strong-on site repul- 
sion U. 

In this paper we describe an attempt to draw conclu- 
sions about the relative importance of U, V and J low- 
energy effective interactions from the numerical arcana of 
cuprate superconductivity, by requiring quantitative con- 
sistency between weak-coupling descriptions of inelastic 
neutron scattering resonance (INSR) and ARPES data in 
Bi2Sr2Cai_,Y,Cu208/Bi2Sr2CaCu208+5 (Bi2212). To 
describe the INSR we derive weak-coupling expressions 
for the (tt, tt) dynamic spin-susceptibility which are exact 
for any U — V — J model. The spin-susceptibility appears 
as a member of a quartet of coupled response functions; 
we find that a four-channel response-function theory is 
required for any non-zero value of J7, V or J. 

We view the use of a weak-coupling theory, in which 
the fermionic excitations are BCS theory Bogoliubov 
quasiparticles, as something which is justified at low tem- 
peratures by experiment; there is nothing manifestly ex- 
otic about the cuprate superconducting state apart from 
its d-wave character and the large energy scales. This 
work is motivated by the expectation that a clearer un- 
derstanding of the character of the low-energy effective 
interactions might hint at its microscopic origin. Starting 



from the assumption that only the U, V and J effective 
interactions are important and using ARPES experimen- 
tal values for the square lattice energy band which crosses 
the Fermi energy, we conclude from ARPES antinodal 
gap values we conclude that in the moderately under- 
doped regime 3J/2 — 2V ~ 250meV. Similarly from 
the ocurrence of the INSR phenomenon we conclude that 
2J + U ^ 350meV. The proximity of these two energy 
scales strongly suggests that the Heisenberg effective in- 
teraction J is dominant. We argue that this finding 
suggests that superconductivity is mediated by short- 
range antiferromagnetic superexchange interactions be- 
tween low-energy quasiparticles which are a remnant of 
the parent antiferromagnetic Mott-insulator and discuss 
some of the challenges which stand in the way of a com- 
pletely satisfactory microscopic theory of superconduc- 
tivity in doped Mott insulators. 



II. PHENOMENOLOGICAL MODEL FOR THE 
SUPERCONDUCTING STATE AT T=0 

A. Effective Hamiltonian 

We consider a low-energy effective Hamiltonian for un- 
derdoped Bi2212, 



(1) 



where e{k) is the band energy and the interaction terms 
are 



Hj — J Si ■ Sj. 

<h3> 



The angle bracket notation < i,j > is used to specify 
that the V and J interactions in our model are restricted 
to near neighbours. We do not view this phenomeno- 
logical Hamiltonian as microscopic, but as what remains 
after incoherent higher energy electronic fluctuations are 
integrated out. In view of the Luttinger theorem, the 
chemical potential /i is nevertheless fixed by the doping 



2 



concentration a; = 1 — y^r: (ct cr )■ The effective inter- 
action parameters [/, V , and J are assummed to be at 
least weakly doping dependent. 



B. Antinodal Gap 

The order parameter for the d-wave superconducting 
state is (ci|Ci+f |) = (— )'^A, where {—Y — +1 for f — ±x 
and —1 for f — ±y, accounting for the d-wave character. 
Applying BCS mean-field theory to Eq.( [T]) leads to a 
gap function A(fc) = VpA (cosfc^ — cosfcj^) and to band 
energies e — s- e' which are modified by interactions. Here 



Vp = 3 J/2 - 2V 



(3) 



has contributions from both of the interactions which 
can induce d-wave superconductivity. Mean-field ground 
state properties are completely determined by e'(fc) and 
Vp. The BCS Hamiltonian yields quasiparticles energies 
±E{k) — ±{e'^{k) + A^(fc))^/^ which are measured in 
ARPES experiments. For e'(fc) we use the experimental 
Bi2212 normal state quasiparticle band structure^. The 
pairing potential Vp of the d-wave superconductor is fixed 
by setting the mean-field maximum gap Ao = 2|V^|A 
equal to the ARPES antinodal gapl°'". Table |T] sum- 
marizes Vp values obtained in this way for several un- 
derdoped Bi2212 samples. For concreteness we focus our 
discussion of numerical consistency between ARPES and 
INSR on the case of doping x = .144, reserving a discus- 
sion of doping dependence to the end of the paper. For 
X = .144, Vp ~ 250meV. One of the central question of 
cuprate superconductivity theory is whether this pairing 
is due primarily to an attractive spin-independent effec- 
tive interaction {V < 0) or primarily due to an antiferro- 
magnetic spin-dependent effective interactions (J > 0). 



C. Competing Orders 

The conclusions reached in this paper depend criti- 
cally on using the same weak-coupling Hamiltonian to 
consistently describe ARPES quasiparticle data and the 
INSR feature in neutron scattering experiments. As we 
explain in more detail below, the emergence of a spin 
resonance well below^^ the particle-hole continuum sig- 
nals an incipient instability in cuprates, almost certainly 
the instability toward the antiferromagnetic state. In a 
weak-coupling generalized random-phase approximation 
(GRPA) theory the energy cost of static antiferromag- 
netic fluctuations is the sum of quasiparticle and in- 
teraction energy contributions. The quasiparticle con- 
tribution K^p is a property of the mean-field state and 
based on ARPES data we can conclude that its value is 
~ 400meV. We find below that =K^p-2J-U and 
conclude from qualitative and quantitative aspects of the 
INSR phenomenon that ivT* <C K^p-, more quantitatively 
a value close to ~ 50meV seems likely. It follows that 



TABLE I: Singlet-pairing potential Vp for several underdoped 
Bi2212 samples. The doping x is extracted from experimental 
Tc data, assuming the parabolic relation proposed by Presland 
et al. in Ref.[L13]. 



X 


n (K) Vp 


(me 


V) fi (meV) 


0.144 


92 


250 


-116.467 


0.126 


85 


256 


-111.358 


0.11 


75 


278 


-105.584 


0.099 


65 


284 


-102.369 



2J + U ^ 350meV. This conclusion is consistent with 
many experiments which hint at a close competition^^ be- 
tween spin and superconducting order in cuprates. To ex- 
plain this assessment more completely, it is necessary to 
describe the weak-coupling theory of spin and supercon- 
ducting fiuctuations in d-wave superconductors in greater 
detail. 



III. WEAK COUPLING INSR THEORY 

Because the interactions in our model Hamiltonian 
are either on-site or nearest-neighbor, fluctuations at the 
high-symmetry wavevector q = Q, where Q is the or- 
dering wavevector of the parent antiferromagnetic Mott 
insulator, can be expressed in terms of a small number 
of coupled channels^. To demonstrate this, we spe- 
cialize to the Sz = +1 projection of the triplet fluc- 
tuation spectrum which is the one relevant to the INS 
measurements. Expanding the Hamiltonian to quadratic 
order around the mean-field state leads to the fluctuation 
Hamiltonian^: 



Hf{t) 



1 
A 



J2 Hk-P) 5{cl^cp.ijs)c\^^-c^, 



p,k,q,(T 



-\-J(k — p) S(ct ct )cr t -\- h.c. 

' \ ^/ \ q—pa pa' ka q — ka 

(4) 



where 

I{k~p)^ [-U + 2V{cos{k^-p.^) + coB{ky-py))l (5) 
and 

J{k - p) = - V[cos{kx -Px)+ cos{ky - py)] . (6) 

The momentum-dependent term which appears in both 
interaction form factors, I{k — p) and J(fc — p), can be 
rewritten as a sum of separable contributions using 

cos{kx-Px)+cos{ky~py) = {si^sp+di^dp;+ss^sSp;+sd^sdf)/2, 

(7) 

where 



(cos kx + cos ky) 
- (smkx + sin/cy) 



d- = (cOS/Ca; - COS ky), 

sd^ = {sinkx — sin ky). 



(8) 
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At wavevector q = Q, the channels involving sine func- 
tions all vanish and we then can identify seven operators, 
some with s- or c?-wave form factors (s^: or d-^), whose 
fluctuations are influenced by interactions: 



is = — ^ V df;S+ 




A, 



-l^Y,^p{Dp + Dp 

p 



(9) 



where Si = ct^c-^^^j^ is a spin-flip operator and 

C0_p^cpi and Djj = cL ^^ct| are pair annihilation and 

creation operators. It follows that ^3 induces d-wave pair 
amplitude and A4 d-wave pair-phase applications. The 
two-particle Greens functions which capture the fluctua- 
tions of these operators are: 



Xa6(Q,w) 



dte^-'6mA,{t),Al{Q)]). (10) 



We focusiS. on the s-wave spin and d-wave pair fields 
(Ai_4), which decouple from the d-wave spin and s-wave 
pair fields (A^-j) and are responsible for the INSR. We 
emphasize that none of the earlier theoretical interpre- 
tations of the INSR feature at q — Q have accounted 
properly for four-channel coupling which appears when 
the weak-coupling calculation is executed correctly. 

The GRPA Greens functions can be obtained by 
solving the equation-of-motion for Xab{Q,^) with the 
quadratic Hamiltonian given in Eq. [4l Since BCS mean- 
field theory with effective interactions describes the low- 
energy charged excitations probed by ARPES, the same 
theory should also describe the low-energy particle-hole 
excitations probed by inelastic neutron scattering. We 
find that 



X 



XqpiQ,^) 



V 



(11) 



where V = diag(-[/ - 2 J, J/2 
the interaction kernel and 



fafb 



-2V,V + J/4, V + J/4) is 



(-l)°+7a/b \ 

E{k) - E{k') uj + Eih) + E(k') ) 

(12) 

is the bare mcan-field-quasiparticle response function. In 
Eq.([T2l) k' = Q~k, the factor {-1)"+'' specifies the sim- 
ple relationship between quasiparticle pair-creation and 
pair-annihilation matrix-elements'' at q = Q, and the 
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FIG. 1: XqpabiQ^^) ^'^^ LO < flo where Qq is the gap 
in the quasiparticle-pair excitation spectrum ai q — Q es- 
tablished by d-wave order. For each channel a, the solid, 
dotted, dashed, and dash-dotted lines represent respectively 

Xqp,ai,Xqp,a2,Xgp,a3' and Xqp,a4- i^o « 70meV for x=0.144 
using the Vp value listed in Table HI) In a weak-coupling 
theory these results depend only on the normal state band 
dispersion and on = 3J/2 — 2V. Channel a — 1 corre- 
sponds to spin, a = 2 to spin with an extended s-wave form 
factor, a = 3 to d-wave pair amplitude and a = 4 to d-wave 
pair phase. 



form factors fa are^: 

/ = ip-{k, k'), ^P^{k, k'),d^l+{k, k'), d^J-{k, k')), 

p^{k,k') = 1 fe' I' i^'' 



V2 



l^(k,k') 



V2 



(13) 



In the GRPA the lo dependence of comes entirely 
from the cj-dependence of Xqp , which depends only on 
the band-structure, on the doping x, and on Vp. Typical 
numerical results for Xqp{Q) are summarized in Fig. [[T]. 

Since the INSR frequency is well below the lower-edge 
of the q — Q particle-hole continuum, it is useful to ex- 
pand Xqp,ab{Q,^) to leading order in lo: 

Xqp,ab{Q,i^) ~ -Ro(a, fe) - Riia,b)uj + O(w^) 

-fafb [1 + (-1)-^+^] 



Ro{a,b) = Y,- 

k 



E{k)+E{k') 

/aA[l^(-l)"+^] 

[E{k)+E{k')Y 



(14) 



The leading coupling between even and odd a operators 
is the Berry-phase coupling which appears at first order 
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in w; the most important^i^ of these is the coupling 
between spin (a = 1) and d- wave-pair phase (a = 4). 
Even-even and odd-odd fluctuations have no Berry phase 
coupHng, but are coupled in the static limit. Xqp has 
a similar low-frequency expansion in which even-even 
and odd-odd fluctuations have relatively little frequency- 
dependence until Lu approaches the particle-hole contin- 
uum closely as seen in Fig.[ [Tj. The even-even and 
odd-odd elements of — = K specify the energy cost 
of the corresponding particle-hole channel fluctuations 
while the even-odd elements, approximately linear in fre- 
quency, specify how the collective fluctuation energy is 
quantized. The even-odd elements satisfy Xqpab ~ ^^ab- 
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IV. MAGNETIC PLASMON 



The INSR energy E''"' solves 
detlr'l =det|x°- 



V\ = 0. 



(15) 



The results for E''''^^ which are implied by Eq[T5]when it is 
assumed that the antifcrromagnetic J interaction is dom- 
inant are summarized in Fig. [[5]; the same values of E^"^^ 
can be obtained by the time-dependent mean-field the- 
ory described in our ealier work^ which did not specialize 
to and take advantage of the simplifications possible at 
q = Q. The agreement of these results with INSR data is 
remarkable. If one interaction is assummed to be domi- 
nant the model parameters are completely determined by 
ARPES data. The Heisenberg effective interaction model 
predicts accurate INSR position values over a broad range 
of doping. 

To achieve a qualitative understanding of this find- 
ing we neglect the (a = 2) extended-s spin-density 
fluctuations which are much stiffer than other fluctua- 
tion modes, as shown in Fig.[[ljb)], and the weak fre- 
quency dependence of the fluctuation energy contribu- 
tions. With these approximations^^ 



E"" 



(16) 

where = Kff - U -2J, K""' ^ K^^ + V + J/4, 
and K'^ — K'^ + V + J/4 are spin, tt amplitude mode, 
and TT phase mode stiffnesses respectively. Since C14 ~ 2 
and the experimental value of E'''^* ^ 40meV is small 
compared to « 400meV and KH « 200meV, it is 
clear that interactions must substantially reduce the val- 
ues of either or JT"^, or both. In the Heisenberg only 
model the interaction J reduces to values which are 
~ 50meV and decrease as expected as the antiferromag- 
nctic state is approached. 

The lines in our three dimensional interaction param- 
eter space which are consistent with both ARPES and 
INSR data are illustrated for several doping values in 
Fig. [ [3] . We note that the intersection of all lines with 
{/ = also yields ^ = to an excellent approximation. 



FIG. 2: INSR energy E'"'" calculated from Eq. [15] (solid 
dots) assumming that with U = V = Q. When only J is non- 
zero, its value (J — ^Vp) is fixed by the ARPES antinodal 
gap. As shown in this figure this assumption predicts the 
correct value for the INSR position. The long-dashed line 
plots in this figure plot the empirical rule = SAIibTc. 

The triangles, white dots, and the black squares show the 
doping dependencies of Oq, and 2Ao calculated from Vp and 
the band structure and the value of K" obtained when these 
are combined with the ARPES determined value J in the 
U = model. 



If we imagine that the U, V and J effective interactions 
have independent origins, it is natural to expect that one 
of the three is likely to be dominant. Our analysis is 
consistent with this expectation and selects the Heisen- 
berg interaction among the three possibilities. This ar- 
gument suggests a remarkably simple single-interaction 
low-energy effective model with strength J ranging from 
« 166 meV for x=0.144 to « 190 meV for x=0.099. 

The INSR position could of course also be accounted 
for by fine-tuning both U and V at fixed Vp^, although 
we have argued that this is less likely. For example, if 
we first assume that d-wave pairing is due entirely to 
attractive spin-independent effective interactions, V — 
-~Vp/2 w — 130meV. This value of V results in a small 
phase stiffness, K'^ ~ 70meV and would require that 
U ~ SOOmeV in order to reduce E^^'^ into the experimen- 
tal range. As we explain later, this parameter set would 
correspond to a weaker correlation scenario in which the 
effective value of U is still fairly large. The INSR mag- 
netic plasmon in this case has approximately equal pair- 
phase and spin character; in the large J scenario on the 
other hand the INSR mode has dominant spin character 
because JC* ^ K'^. None of these arguments are suffi- 
cently quantitative to favor the V — model choice over 
the V = —J/4 choice, commonly used in i — J model 
calculations and motivated by the theory of the superex- 
change mechanism. Compared to the J-only model, this 
choice shifts J slightly, from J = 2K,/3 to Vp/2, resulting 
in slightly larger and smaller without shifting the 
INSR position significantly. 
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FIG. 3: Lines in {U, J) parameter space which reproduce both 
the experimental antinodal gaps Aq and the INSR energy 
for a series of doping values. Note that V = 3 J/4 — Vp/2 
where the Vp values are listed in Table U For all doping 
values the intersection of these lines with U — also yields 
V = 0. We argue in the text that this is strongly suggestive 
that the Heisenberg J interaction is dominant. This figure 
was constructed using the 4 by 4 expression given in Eq. [15] 
for the resonance frequency. 



We close this section by commenting briefly on the 
role of inter-channel Berry phase coupling, which has of- 
ten been neglected in theoretical analyses, in determin- 
ing the INSR position and character. When only the 
spin-channel is included the equation for the resonance 
frequency is K^^E"^^^) = 0. Because of the weak energy 
dependence of K^^ below the particle-hole continuum we 
see that when coupling is ignored, [/ -I- 2 J has to be ad- 
justed to more than 90% of K^^ to explain the resonance 
position, placing the system even closer to the antifer- 
romagnetic state instability point. For a given value of 
the interaction stregth, mode coupling substantially low- 
ers the resonance frequency. Mode coupling is therefore 
important in explaining the experimental relationship be- 
tween the value of the resonance frequency and the prox- 
imity of the antiferromagnetic state. 



V. DISCUSSION 

This work starts from the observation that in the low- 
temperature limit cuprates act like ordinary supercon- 
ductors with elementary quasiparticle excitations that 
are describable by BCS theory. ARPES experiments find 
a Fermi surface with area proportional to electron dop- 
ing 1 — a;, rather than hole-doping x and a Fermi ve- 
locity which shows no sign^^ of declining with x. The 
most remarkable feature of cuprate superconductors at 
low-temperatures is the sharp collective excitation which 
appears at wavevectors near Q = (tt, tt) in neutron- 
studies of quantum spin-fluctuations. The apparently 



banal character of the low-energy excitations suggests 
that the ground state of these superconductors can be de- 
scribed using a weak-coupling theory in which the quasi- 
particles seen in ARPES experiments interact in a way 
which leads to both d-wave superconductivity and the 
inelastic neutron scattering mode. Given this starting 
point, it is natural to ask if quantitative comparison with 
experiments gives any indication as to the character of 
the presumably renormalized interactions which appear 
in this low-energy effective theory. In making this as- 
sessment, it is important to acknowledge that the prop- 
erties of cuprate superconductors are not entirely uni- 
versal, and not all experimental data is available on any 
one material. Nevertheless we judge that well established 
broad trends allow generic conclusions. The numbers 
mentioned below are intended to apply most closely to 
the Bi2212 cuprate family. 

In order to make progress toward identifing the renor- 
malized interactions it seems to be necessary to start 
by limiting the forms which can reasonably be con- 
sidered. The fact of singlet d-wave superconductivity 
on the cuprate square lattice appears to require either 
spin-independent attractive interactions (possibly due to 
phonons) or spin-dependent antiferromagnetic Heisen- 
berg interactions (possibly due to short-range antiferro- 
magnetic order) between electrons on neighbouring sites. 
The proximity of a Mott insulator state for hole-doping 
X = would seem to require allowing for the possibility of 
a strongly repulsive on-site Hubbard-like interaction U . 
We have therefore considered how the properties of the 
system depend on all three parameters C/, V , and J. The 
simplest experimental constraint on this parameter set is 
imposed by the magnitude of the d-wave antinodal gap 
which, when combined with the band-structure known 
from ARPES, constrains Vp = 3 J/2 — 2V to a value 
^ 250meV which increases moderately with decreasing 
X as summarized in Table I. 

An understanding of the implications of the strong 
INSR at Q = (tt, tt) requires a more subtle analysis. 
In a weak-coupling theory we can determine collective 
excitation properties by examining the influence of in- 
teractions on single particle- hole excitations. As ex- 
plained in Section III, for the case Q = (tt, tt) exci- 
tations of a d-wave superconductor on a square lattice 
the spin-response function probed by neutrons is coupled 
to pair-magnitude, pair-phase, and extended-s spin re- 
sponse functions. A correct analysis of weak-coupling col- 
lective fluctuations therefore requires that all four modes 
be treated simultaneously, rather than singling out spin- 
response^ i^°'^^i^^'^'^i^^'^^i^^ as has been common in RPA 
theories or including only flrst few modes^iii^liS^. We 
emphasize that a theory of the spin-response function 
requires consideration of the four-coupled modes when- 
ever any of the three interactions we consider is non- 
zero. In this theory the energy cost of antiferromag- 
netic spin fluctuations is reduced from Kff ~ 400meV 
to = Kf^ — [/ — 2 J by interactions. {Kff is purely 
a property of the BCS mean-fleld theory and is there- 
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fore specified by the ARPES bands and Vp.) In order for 
the INSR to appear, as it does in experiment, below the 
particle-hole continuum must be reduced to a small 
value by interactions, i.e. the system must be close to 
an antiferromagnetic instability. (A theory which did not 
accout for coupling between spin and other modes would 
require that K'^ be reduced to an even smaller value in 
order to account for the INSR data.) On this basis we 
conclude that U + 2J ^ 350meV. 

Both ARPES and INSR experiments are there- 
fore consistent with a J-dominated renormalized in- 
teraction model. The strength for J is comparable 
to the strength of the superexchange interactions in 
the undoped antiferromagnetic insulator state of the 
cuprates. The effective interactions that are consis- 
tent with experiment are therefore similar to the ef- 
fective interactions which appear— in < — J models, in- 
tended to describe quasiparticlcs which live nearly en- 
tirely in the Gutzwiller projected^^iH-lliHiS no-double- 
occupancy Hubbard-model subspace. The absence of 
a large Hubbard replusion U in the effective interac- 
tion is also expected for quasiparticles which avoid this 
strong interaction by residing dominantly in the no- 
double-occupancy subspace. Any other choice of inter- 
action model requires improbable fine-tuning to achieve 
consistency with the main experimental facts. 

The consistency of the J-only interaction model with 
ARPES and INSR data partially affirms one of the main- 
pictures^ of cuprate superconductivity, namely that it 
emerges from the antiferromagnetic interactions in the 
parent x = Mott insulator state when holes in the 
lower Hubbard band become itinerant at sufficiently 
large doping. There are, however, at this point ma- 
jor challenges which still stand in the way of a fully 
consistent theoretical picture to accompany this sce- 
nario. Mean-field-theory slave-boson and other related 
formulations^ii^ of the t — J model, for example, tend 
to lead to quasiparticle band widths which decreases 
wtih hole-doping and to superfiuid density which is 
strictly proportional to hole-doping x. These predic- 
tions are, however, clearly at odds with ARPES"'^^ and 
penetration-depth2ii2^i^ temperature-dependence obser- 
vations, both of which indicate that the Fermi veloc- 
ity has a weak doping dependence and that the pene- 
tration depth vanishes at a rather large value of hole- 
doping. A second theoretical approach motivated by 
the t — J model proceeds from the ansatz that the 
superconducting state can be described by Gutzwiller- 
projected BCS paired statea2Si^>^i^i2ii^. Variational 
Monte Carlo calculations^ have demonstrated that the 
Gutswiller-projected BCS wavefunction reproduces a va- 
riety of quantities observed in experiments. Furthermore, 
using the Gutzwiller approximation^,the effect of the 
Gutzwiller projector can be replaced approximately by 
doping-dependent band energy and interaction renormal- 
ization factors, leading to mean-field equations similar to 
those of slave-boson mean-field theories and of the same 
form as those of the weak-coupling BCS theory. This 



simple approximation can reproduce the results of the 
variational Monte Carlo studies to a remarkable degree. 
As a result, the plain vanilla theory helps provides a the- 
oretical explanation for weak-coupling low-energy behav- 
ior in the superconducting state with effective interaction 
parameters that arise in a complex way from the corre- 
lated fluctuations of higher energy degrees of freedom. 
From this point of view, our main finding in this paper 
is that the same interactions are also quantitatively con- 
sistent with the INSR feature. In addition, as we have 
shown— previously, the pairing wavevector dependence of 
these low-energy antiferromagnetic fluctuations, provides 
a natural explanation for the fact that the superfiuid den- 
sity goes to zero at a finite value of hole-doping. 

One very interesting property of cuprate superconduc- 
tivities is that the Fermi velocity remains nearly constant 
throughout the underdoped regime, even as^ the su- 
perfiuid density and the low-frequency infrared spectral- 
weight weight decline, presumably reflecting in part a 
decrease in the quasiparticle renormalization factor Z. 
Although this property is reproduced by Gutswiller pro- 
jected variational Monte-Carlo calculations^, it seems 
surprising from a perturbation theory point of view that 
Z should decline to a small value without a corresponding 
decrease in Fermi velocity; large values for the energy- 
dependence of the self-energy can easily arise from small 
energy denominators associated with resonances near the 
Fermi energy, but it is not so obvious how correspond- 
ingly large values for the wavevector-dependence of the 
self-energy could arise. We suspect that the quasiparticle 
renormalization factor in the cuprates is not as small as 
suggested by the Gutzwiller-projected BCS state ansatz. 

Another very interesting property of cuprate supercon- 
ductores is the marginal Fermi liquid behavior which of- 
ten occurs when superconductivity is supressed by tem- 
perature (in the overdoped regime) or by a magnetic field. 
Marginal Fermi liquid behavior suggests that quasiparti- 
cles are strongly coupled to very-low-energy bosonic ex- 
citations, or that the system would be very close to a 
quantum phase transition if superconductivity was sup- 
pressed. It seems likely to be more than a coincidence 
that in the low-energy effective model we have extracted 
from experiment, the cuprates would be very close to an- 
tiferromagnetism throughout much of the superconduct- 
ing dome if they did not form a superconducting state. 
This statement is quantified in the appendix in which we 
discuss quantitative aspects of the competition between 
superconductivity and antiferromagnetism. Decreasing 
the gap Aq ^ decreases Kff by approximately 50 meV, 
just the amount required to reach the antiferromagnetic 
instability. Because superconductivity supresses antifer- 
romagnetism, pair-condensation must raise a part of the 
electron-electron interaction energy. This property may 
be part of the explanation for the apparent increas a'^^i^° 
in kinetic energy in the superconducting state. 

In conclusion, we have performed a weak-coupling 
analysis of ARPES and INSR experiments in Bi2212, 
in an effort to indentify an effective interaction model 
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given approximately by: 
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FIG. 4: Numerical evaluation of —Xqpiii^ = Q,uj = 0) as a 
function of Aq for the doping x = 0.144 using Eq. [121 



which is consistent with both experiments. We find that 
the doping dependences of the superconducting gap, and 
the INSR energy E'^'^'^ can be consistently explained by 
a model with near-neighbor Heisenberg interactions with 
a strength that is consistent with superexchange interac- 
tions. This result suggests that strong short-range repul- 
sion and incoherent remnants of the antiferromagnetic 
insulating parent compound are key to high-temperature 
superconductivity. 
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APPENDIX A: COMPETITION BETWEEN 
d-WAVE SUPERCONDUCTIVITY AND 
ANTIFERROMAGNETISM 



In this appendix we explain in more detail how d-wave 
superconductivity and antiferromagnetism compete in a 
weak coupling GRPA theory. As pointed out in an earlier 
paper™, this competition is responsible for a correlation 
induced suppression of the superfluid density as antifer- 
romagnetism is approached. 

The quantum zero point energy associated with the 
collective fluctuations probed in INSR experiments is 



1 ' 

E'^P ^ -Y^ El"'' = 



Ki 

k.adia 



Ki 

k.adia 



c 



(Al) 



14. fc 



where the prime refers to the sum over wavevectors near 
k = Q for which the Berry curvature C14 is large. 
Ki , . , Ki are adiabatic limit stiffnesses which equal 

k,adia k,adia 

-X-^n(fc, w = 0) - [/ - 2J and -x^p.^ik^uj = 0) + V + 
J/4 respectively. The correlation contribution to the su- 
perfluid density is then given bj*^: 



1 dE^'P 
A 9P2 



I 



J^res 
k 



1 

k k,adia 



k.adia 



Qp2 



(A2) 



where P is the norm of the pairing wavevector of the 
BCS wavefunction. (We explain below why the P de- 
pendence of K'^ is much stronger than the P dependence 
oi K't'.) Since both El" 

k 



and Ki , are positive quan- 

k.adia 



titles, dKi ^ IdP^ 



■'^ determines the sign of p" 
the chain rule leads to 

dKf. 



k.adia 



ap2 



ap2 dA, 







Using 



(A3) 



Because we know that d-wave superconductivity is weak- 
ened by pair-breaking effect of finite P Doppler shifts on 
nodal quasiparticles we can conclude that dAo/dP'^ < 0. 
A positive value for dKi Id An characterizes com- 

k^adia 

petition between antiferromagnetism and superconduc- 
tivity. When the two orders compete, the energy cost 
of antiferromagnetic fluctuations is increased when su- 
perconductivity strengthens. In other words, will 
increase as Aq increases, giving dKf /OAq > in 

this case. To confirm the competing nature of these in- 
teractions, we have performed a numerical calculation 
of Aq dependence of K'^. Since K'^ — —Xqp^ii{k,uj = 
0) — [/— 2 J, the Aq dependence of K'^ comes entirely from 
~Xgp,iiik,uJ = 0). Fig. [H] plots -x-p\,ik = Q,u; = 0) 
as a function of Aq for doping x = 0.144 as an exam- 
ple. —Xqp iiik = Q,uJ = 0) clearly increases monotoni- 
cally with Ag, consistent with the competing order pic- 
ture discussed above. It is important to realize that the 
relative change in is larger than the relative chance 
in —Xqp iiik = = 0) by apprxoximately a factor of 
ten because of the interaction contributions to the inverse 
response functions. We have shown separately* that the 
size of this effect combined with the reduction in Aq at 
finite P can explain the large reduction in superfluid den- 
sities compared with mean-field theory values in under- 
doped cuprates. 

Our demonstration that superconductivity and anti- 
ferromagnetism complete in the effective weak-coupling 
picture of cuprate superconductors contradicts the the 
conclusion of Zhang et alM who state that weak-coupling 
theory predicts an enhancement of antiferromagnetism 
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due to d-wave superconductivity. Our calculation shows 
that their conclusion is not valid when the the nor- 
mal state Fermi surface is similar to what is observed 
in Bi2212 systems. As pointed out by Tchernyshyov 
et al., the response functions x*^ at q* = Q are dom- 
inated by the hot spots located near (tt, 0). If this is 
true, -Xqp^iiik,uj = 0) ~ 2Aantinodai = 2Ao. This 
crude argument is consistent with dKt , /dAn > 0, 

k.adia' 



and very crudely consistent with Fig. [[4]. The character 
of the competition between superconductivity and anti- 
ferromagnetism is dependent on the normal state band 
structure. This dependence is likely responsible for some 
of the differences between hole-doped and electron-doped 
cuprates. 
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